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We investigate constraints on Lorentz invariance violation in the neutrino sector from a joint 
analysis of big bang nucleosynthesis and the cosmic microwave background. The effect of Lorentz 
invariance violation during the epoch of big bang nucleosynthesis changes the predicted helium- 
4 abundance, which influences the power spectrum of the cosmic microwave background at the 
recombination epoch. In combination with the latest measurement of the primordial helium-4 
abundance, the Planck 2015 data of the cosmic microwave background anisotropies give a strong 
constraint on the deformation parameter since adding the primordial helium measurement breaks 
the degeneracy between the deformation parameter and the physical dark matter density. 


PACS numbers: 98.80.Es 


1. INTRODUCTION 


Neutrino oscillations have shown that there are small but nonzero mass squared differences between three neutrino 
mass eigenstates (see Ref. and the reference therein), which imply the existence of physics beyond the standard 
model of particle physics. However, neutrino oscillation experiments cannot tell us the overall mass scale of neutrinos. 
Fortunately, cosmology provides a promising way to determine or constrain the total mass of neutrinos by the gravi- 
tational effects of massive neutrinos which can significantly change the CMB power spectrum due to lensing [2H5] and 
the formation of large-scale structures [6] [7], therefore, alter the cosmic microwave background (CMB) anisotropies 
and the large-scale structure distribution of matter (see [8] for a review). Recent Planck 2015 data combined with 
some low-redshift data give an upper limit on the total mass of neutrinos at 95% confidence level, Nm, < 0.19 eV [5]. 

Another possible signal of new physics is Lorentz invariance violation in the neutrino sector. The observed neutrino 
oscillations may originate from a combination of effects involving neutrino masses and Lorentz invariance violation [9}+ 
[14]. Lorentz invariance is a fundamental symmetry in the standard model of particle physics. Although present 
experiments have confirmed Lorentz invariance to a good precision , it may be broken in the early Universe when 
energies approach the Planck scale. The standard model itself is believed to be a low-energy effective theory of an 
underlying unified theory. Lorentz invariance violation has been explored in quantum gravity [I6], loop quantum 
gravity , hon-commutative field theory , and doubly special relativity theory [I9]. Various searches for Lorentz 
invariance violation have been performed with a wide range of systems Bil. 

The cosmological consequences of neutrinos with the Coleman-Glashow type [I0] dispersion relation have been 
studied in Ref. [22]. As can be expected, the Lorentz-violating term affects not only the evolution of the cosmological 
background but also the behavior of the neutrino perturbations. The former changes the expansion rate prior to 
and during the epoch of photon-baryon decoupling, which alters heights of the first and second peaks of the CMB 
temperature power spectrum, while the latter alters the shape of the CMB power spectrum by changing neutrino 
propagation. Since these two effects can be distinguished from a change in the total mass of neutrinos or in the 
effective number of extra relativistic species, CMB data have been proposed as a probe of Lorentz invariance violation 
in the neutrino sector. 

Moreover, since the Lorentz-violating term influences the abundances of the light elements by altering the energy 
density of the Universe and weak reaction rates prior to and during the Big bang nucleosynthesis (BBN) epoch, 
BBN provides a promising probe of Lorentz invariance violation in the neutrino sector in the early Universe |23|. 
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In particular the BBN-predicted abundance of helium-4 is very sensitive to the Lorentz-violating term. It is well 
known that the primordial helium-4 abundance plays a non-negligible role at the recombination epoch because it 
has an influence on the number density of free electron. Therefore, Lorentz invariance violation affects the CMB 
power spectrum by changing the BBN-predicted abundance of helium-4. A reasonable way to test Lorentz invariance 
violation in the neutrino sector with CMB data is to take account of the BBN-predicted helium-4 abundance as a 
prior rather than fixing the helium-4 abundance. 

In this paper we investigate constraints on Lorentz invariance violation in the neutrino sector from a joint analysis 
of BBN and CMB. In combination with the latest measurement of the primordial helium-4 abundance, the Planck 
2015 data of the CMB anisotropies give a strong constraint on the deformation parameter. Adding the primordial 
helium measurement can break effectively the degeneracy between the deformation parameter and the physical dark 
matter density. 

This paper is organized as follows. In Sect. [2] we parameterize Lorentz invariance violation in the neutrino sector. 
In Sect. |3| we derive the weak reaction rates in the Lorentz-violating extension of the standard model and calcu- 
late the BBN prediction of the helium-4 abundance. In Sect. |4} we derive the Boltzmann equation for neutrinos in 
the synchronous gauge and calculate the CMB power spectrum. In Sect. |5| we place constraints on the deforma- 
tion parameter using the Planck 2015 data in combination with the latest measurement of the primordial helium-4 
abundance. Sect. |6]is devoted to conclusions. 


2. DEFORMED DISPERSION RELATION 


At a phenomenological level, the deformed dispersion relation can be characterized by a power series of momentum. 
In this paper, we consider a simple instance of the dispersion relation for neutrinos constructed in [10], which can be 
parameterized as follows: 

E? =m? +p? + Ep, (1) 
where F is the neutrino energy, m the neutrino mass, p = (p'p;) 12 the magnitude of the 3-momentum, and € the 
deformation parameter characterizing the Lorentz symmetry violation. 

In the spatially flat Friedmann-Lemaintre-Robertson-Walker metric with the scale factor a, the number density ny, 
energy density p, and pressure P, for neutrinos with are given by 
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The phase space distribution for neutrinos is fo(q) = gs [1 + exp (€/To)] 7}, where gs = 2 is the number of spin 
degrees of freedom, To the neutrino temperature today, q the magnitude of the comoving 3-momentum, and € = 
/m?a? + (1 + £) q? the comoving energy. 


3. BBN PREDICTION 


The abundances of the light elements produced during the BBN epoch depend on the competition between the 
nuclear and weak reaction rates and the expansion rate of the Universe. In the standard cosmological scenario and in 
the framework of the electroweak standard model, the dynamics of this phase is controlled by only one free parameter, 
the baryon to photon number density. 

We consider N species of nuclides whose abundances X; are the number densities n; normalized with respect to 
the baryon number density ng, 
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Their evolutions are ruled by the following Boltzmann equations, 
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where I’ denotes the reaction rate and N; is the number of nuclide 7 involved in the reaction. 
In the case of Lorentz invariance violation in the neutrino sector, the changes of the reaction rates of the following 
weak reactions have to be taken into account: 


nt+ve @ pte, (7) 
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which determine the neutron-to-proton ratio when the baryons become uncoupled from the leptons. The abundances 
of the light elements depend on the neutron-to-proton ratio at the onset of BBN. As an example, let us compute the 
reaction rate of the process in the Lorentz-violating extension of the standard model [20] [24]. The reaction rate 
is 
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where Cy and Cy are the vector and axial coupling of the nucleon, P) is the standard reaction rate per incident 
nucleon derived in Ref. [25]. The first prefactor on the right-hand side of Eq. arises from the neutrino propagator 
and the evW coupling in the Lorentz-violating extension of standard model. The second prefactor comes from the 
statistical distribution for neutrinos. As for the other processes in (7)-(9). their reaction rates can be simply derived 
by properly changing the statistical factors and the delta function determined by the energy conservation for each 
reaction. Therefore, the corrections to the conversion rate of neutron into proton and its inverse rate are the same as 
in (20). 

In order to calculate the abundance of the light elements produced during the BBN epoch, we modified the publicly 
available PArthENoPE code to appropriately incorporate the Lorentz-violating term in the neutrino sector. As 
shown in [23], the BBN-predicted abundance of helium-4 is sensitive to the deformed parameter € due to the fact that 
it is mainly determined by the neutron-to-proton ratio, which is related to the weak reaction rate, neutrino number 
density and expansion rate of the Universe. In our analysis we shall focus on the helium-4 abundance that influences 
the CMB power spectrum. 

Although in our analysis we use an interpolation method to obtain the value of Y, on a grid of points in (Qsh?, €) 
space (see Sect.|5), it is possible to describe the dependencies of the BBN-predicted helium-4 abundances upon the two 
parameters by simple, linear fits which, over their ranges of applicability, are accurate to a few percent. Following [27], 
we get the following linear fit for Y, versus Qh? and £: 


Yo) = 0.2334 + 0.520 4h? + 0.258€, (11) 


over our adopted range in Q,h? and £, as shown in Fig.|1| We note that this fit is a very good approximation over the 
adopted parameter ranges, which is useful in facilitating studies of the viability of various options for non-standard 
physics and cosmology, prior to undertaking detailed BBN calculations 

Since the neutrino species share the temperature with the photons prior to neutrino decoupling, the prediction for 
the effective number of neutrino species is 


Nes = 3.046(1 + €)~°/?. (12) 


The total relativistic energy density p, in neutrinos is given in terms of the photon density p, by 
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4. CMB POWER SPECTRUM 


The Lorentz-violating term affects not only the evolution of the cosmological background but also the behavior of 
the neutrino perturbations. Following Ref. [28] we shall derive the Boltzmann equation of the distribution function 
for neutrinos with (1) to linear order in the perturbations. In the synchronous gauge, the perturbed energy density 
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FIG. 1: BBN-predicted curves for Yp in Qh? (red) and € (blue). From bottom to top, the red curves correspond to € = —0.1, 
0, —0.1, while the blue curves correspond to Nh? = 0.01, 0.02, 0.03. The dotted curves are the BBN-predicted results, while 
the solid curves are our fits. 


py, pressure 0P,, energy flux 0,, and shear stress o, are, respectively, given by 


ion = i f EEAO, (14) 
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where W; are the expansion coefficients of the perturbed neutrino distribution function expanded in a series of Legendre 
polynomials 
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where 7 is the conformal time. These perturbations WV; satisfy the following Boltzmann equations: 
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in the synchronous gauge, where dots denote derivatives with respect to the conformal time, h and ņ are the two scalar 
modes of the metric perturbations in the Fourier space k. To avoid the reflections from high-@ equations by simply 
set Wp = 0 for l > Lmax, such a Boltzmann hierarchy is effectively truncated by adopting the following scheme [28]: 
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Making a transformation q > y1 + &q, we note that the number density, energy density and pressure of neutrinos 
are proportional to (1 + €)~%/?, which means that increasing € decreases the number density, energy density and 
pressure of neutrinos. Moreover, there is a factor \/1 + £ in the second term of the Boltzmann equations, which comes 
from neutrino propagation. 


For the massless neutrinos, the Boltzmann equations are simplified by setting € = /1 + £q. The number of variables 
can be reduced by integrating out the q dependence in the neutrino distribution function so that 
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Using the orthonormality and recursion relation of the Legendre polynomials, from the Boltzmann equations we 
derive the evolution equations of the perturbed energy density 6, = 6p,/p , energy flux 6,, shear stress o, independent 
of the momentum q, 
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The truncation scheme for massless neutrinos is 
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In the synchronous gauge the adiabatic initial conditions for the metric perturbations and massless neutrinos in the 
cold dark matter frame are 
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where C is a dimensionless constant determined by the amplitude of the fluctuations from inflation, Ry = p,/(p,+pv), 
and w = QmHo/ yQ + Qu, which corresponds to the matter contribution to the total energy density of the Universe. 
In the radiation-dominated era the massive neutrinos are relativistic. The initial condition of Yo for massive neutrinos 
is related to 6,. Then, using Eqs. and ignoring the mass terms in the differential equations of high-order moment, 
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FIG. 2: The matter power spectrum for € = —0.1 (dashed), 0 (dash-dotted), 0.1 (solid). 


we get the initial conditions of Y; for massive neutrinos: 
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In order to compute the theoretical CMB power spectrum, we modified the Boltzmann CAMB code to ap- 
propriately incorporate the Lorentz-violating term in the neutrino sector. As pointed out in [22], the effects of the 
Lorentz-violating term on the CMB power spectrum are distinguished from a change in the total mass of neutrinos 
or in the effective number of extra relativistic species [B0H34]. Therefore, the measurements of the CMB anisotropies 
provide a cosmological probe of Lorentz invariance violation in the neutrino sector. Moreover, we calculate the matter 
power spectrum today for € = —0.1,0,0.1. From Fig. [2] we find that increasing € can enhance the matter power spec- 
trum at small scales. It is well known that variances in the helium-4 abundance modify the density of free electrons 
between helium and hydrogen recombination and therefore influence the CMB power spectrum. However, the effects 
of the Lorentz-violating term on the BBN-predicted helium-4 abundance were not considered in [22]. In the next 
section, we shall put constraints on the deformation parameter from a joint analysis of BBN and CMB. 


5. JOINT ANALYSIS FROM BBN AND CMB 


In our analysis we use a modified version of the publicly available CosmoMC package to explore the parameter space 
by means of Monte Carlo Markov chains technique [85]. Since the deformation parameter € is nearly uncorrelated 
with the total mass of neutrinos when CMB data are used to detect the signature of Lorentz invariance violence [22], 
we only consider massless neutrinos with vanishing chemical potentials in our analysis. Our cosmological model is 
the spatially flat ACDM plus three types of massless neutrino with the deformed dispersion relation (1), which can 
be described by the following seven parameters: 


{Qyh?, Qh’, OMC, Tre; Ns, As, E} , 


where h is the dimensionless Hubble parameter defined by Hp = 100h km s~! Mpc™!, Qh? and 0h? are the physical 
baryon and dark matter densities relative to the critical density, @yc is an approximation to the ratio of the sound 
horizon to the angular diameter distance at the photon decoupling, Tre is the reionization optical depth, ns and A, 
are the spectral index and amplitude of the primordial curvature perturbations at the pivot scale ky = 0.05 Mpc7!. 


TABLE I: Mean values and marginalized 68% confidence level for the deformation parameter and other cosmological parameters, 
derived from the Planck data in combination with baryon acoustic oscillation data, the JLA sample of Type Ia supernovae and 
the measurement of the helium-4 mass fraction. For comparison with BBN consistency, fixing Yp = 0.24 is also considered. 


Model Planck2015+JLA+BAO Planck2015+JLA+BAO+Y, 
Yp = 0.24 (fixed) | BBN consistency BBN consistency 
wh? 0.02127 + 0.00046 | 0.02212 + 0.00025 0.02227 + 0.00021 
Qh? 0.1124 + 0.0061 | 0.1163 + 0.0040 0.1197 + 0.0018 
1008mc 1.04181 + 0.00098 1.04173 + 0.00101 1.04090 + 0.00050 
Tre 0.071 + 0.018 0.076 + 0.018 0.077 + 0.018 
ns 0.9583 + 0.0069 | 0.9679 + 0.0049 0.9659 + 0.0044 
lIn(10*° As) 3.064 0.038 3.084 + 0.035 3.087 + 0.035 
Ho(km s~*Mpc')| 65.92 + 1.76 67.25 + 1.00 68.01 + 0.63 
os 0.840 + 0.014 0.846 + 0.016 0.841 + 0.015 
Age(Gyr) 14.191 + 0.35 13.937 + 0.19 13.767 + 0.068 
Qn 0.691 + 0.008 0.694 + 0.007 0.693 + 0.007 
ie 0.309 + 0.008 0.306 + 0.007 0.307 + 0.007 
£ 0.101 + 0.091 0.042 + 0.049 —0.002 + 0.014 


We firstly run the PArthENoPE BBN code with two free input parameters, 0,h? and £, to precalculate the 
primordial helium-4 abundance on a grid of points in (Q,h?, £) space. For each couple of parameter values the helium- 
4 abundances are then obtained from the grid by two-dimensional cubic spline interpolation, which is adopted to 
calculate the CMB power spectrum. This method is called “BBN consistency”. 

We use the recently released Planck 2015 likelihood code and data, including the Planck low-é likelihood at mul- 
tipoles 2 < £ < 29 and Planck high-£ likelihood at multipoles £ > 30 based on pseudo-C, estimators [5]. The former 
uses foreground-cleaned LFI 70 GHz polarization maps together with the temperature map obtained from the Planck 
30 to 353 GHz channels by the Commander component separation algorithm over 94% of the sky. The latter uses 
100, 143, and 217 GHz half-mission cross-power spectra, avoiding the galactic plane as well as the brightest point 
sources and the regions where the CO emission is the strongest. “Planck 2015” denotes the combination of the low-é 
temperature-polarization likelihood and the high- temperature likelihood. 

We use the “Joint Light-curve Analysis” (JLA) sample of type Ia supernovae, which is constructed from the 
SNLS and SDSS supernova data, together with several samples of low-redshift supernovae [36|. Baryon acoustic 
oscillation (BAO) measurements are another important astrophysical data set, which are powerful to break parameter 
degeneracies from CMB measurements. We use BAO measurements of Dy/Trarag from the 6dFGS at zef = 0.106, 
SDSS Main Galaxy Sample at zef = 0.15, BOSS LOWZ at zef = 0.32, and BOSS CMASS at zeq = 0.57 [5]. Here 
Dy is the effective distance measure for angular diameter distance, farag is the comoving sound horizon at the end 
of the baryon drag epoch and Zeg is the effective redshift. In our analysis we also use the helium-4 mass fraction 
measurement of Y, = 0.2449+0.0040, derived from helium and hydrogen emission lines from metal-poor extragalactic 
HII regions and from a regression to zero metallicity [87]. 

In Table [I] we list our results. With the combined data of Planck 2015+JLA+BAO, the deformation parameter 
is estimated to be € = 0.042 + 0.049 if the helium-4 abundance is set from the BBN prediction. For comparison, 
if the helium-4 abundance is not varied independently of other parameters, with the same data set we find € = 
0.101 + 0.091, which indicates that the constraint of € is not improved compared to the result derived in from 
the 7-year WMAP in combination with lower-redshift measurements of the expansion rate. The reason is that 
the estimated values of cosmological parameters are biased by fixing Yp = 0.24. The gray and red regions in Fig. 
show 68% and 95% contours in the €-0,h? plane setting Yp = 0.24, derived from Planck 2015+JLA+BAO and from 
WMAP7+H,+BAO [22], respectively. We see that the Planck 2015 data marginally favor a positive £ compared 
to the 7-year WMAP data. The addition of the latest measurement of the helium-4 abundance leads to a strong 
constraint on the deformation parameter € = —0.002 + 0.014. As shown in Figure [3] adding the Y, data breaks the 
degeneracy between the deformation parameter and the physical dark matter density. This result shows that no signal 
of Lorentz invariant violation is detected by the joint analysis of CMB and BBN. Figure [4] shows the constraints in 
the € — Qh? plane, £ — Hp plane, € — n, plane and € — A, plane from Planck 2015+JLA+BAO. 


6. CONCLUSIONS AND DISCUSSIONS 


We have studied the cosmological constraints on Lorentz invariance violation in the neutrino sector by a joint analysis 
of CMB and BBN. Instead of fixing the value of the helium-4 abundance, we have applied the BBN prediction of 
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FIG. 3: Marginalized two-dimensional contours in the € — Qh? plane with 68% and 95% confidence level in the ex- 
tended ACDM model with BBN consistency, derived from the data of Planck 2015+JLA+BAO (green) and Planck 
2015+ JLA+BAO+Y, (blue). The red and gray contours correspond to Yp = 0.24, derived from Planck 2015+JLA+BAO 
and from WMAP7+H,+BAO [23], respectively. 


the helium-4 abundance determined by € and Q,h? to calculate the CMB power spectra. Using the Planck 2015 
data in combination with the JLA sample of type Ia supernovae and the BAO feature, we have put a constraint on 
the deformation parameter. Adding the measurement of the helium-4 abundance breaks the degeneracy between the 
deformation parameter and the physical dark matter density and so improves the constraint by a factor 3. 

As found in [22], the deformation parameter is nearly uncorrelated with the total mass of neutrinos Um, when 
WMAP data are used to constrain €. Planck 2015 data is much more sensitive to the total mass of neutrinos than 
WMAP data. Therefore, we consider the effects of the neutrino masses on our results. Figure [5] shows 68 and 95% 
contours in the € — im, plane for Planck 2015 data combined with lower-redshift measurements. We see that € is 
nearly uncorrelated with Ym, from Fig. 

Compared to cosmological bounds on the Lorentz-violating coefficient, observations of high-energy astrophysical 
neutrinos give stronger constraints. As listed in Table XIII of Ref. [13], the coefficient is constrained down to O(10~°) 
from the time-of-flight measurements, under the assumption that neutrino oscillations are negligible. Cohen and 
Glashow have argued that the observation of neutrinos with energies in excess of 100 TeV and a baseline of at least 
500 km allows us to deduce that the Lorentz-violating parameter is less than O(10~ 1") [39]. The present work provides 
a new way to probe the signal of Lorentz invariance violation in the early Universe, which can in principle be used to 
constrain £ in the sterile neutrino sector [40/41]. Although present cosmological data are too weak to yield competitive 
constraints, future measurements of CMB and BBN offer prospects for placing stringent constraints. 
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